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ABSTRACT 

WASP-80b is a missing link in the study of exo-atmospheres. It falls between the warm 
Neptunes and the hot Jupiters and is amenable for characterisation, thanks to its host 
star’s properties. We observed the planet through transit and during occultation with 
Warm Spitzer. Combining our mid-infrared transits with optical time series, we find 
that the planet presents a transmission spectrum indistinguishable from a horizontal 
line. In emission, WASP-80b is the intrinsically faintest planet whose dayside flux 
has been detected in both the 3.6 and 4.5 pm Spitzer channels. The depths of the 
occultations reveal that WASP-80b is as bright and as red as a T4 dwarf, but that 
its temperature is cooler. If planets go through the equivalent of an L-T transition, 
our results would imply this happens at cooler temperatures than for brown dwarfs. 
Placing WASP-80b’s dayside into a colour-magnitude diagram, it falls exactly at the 
junction between a blackbody model and the T-dwarf sequence; we cannot discern 
which of those two interpretations is the more likely. WASP-80b’s flux density is 
as low as GJ436b at 3.6 /rm; the planet’s dayside is also fainter, but bluer than 
HD 189733Ab’s nightside (in the [3.6] and [4.5] b'pjfeer bands). Flux measurements on 
other planets with similar equilibrium temperatures are required to establish whether 
irradiated gas giants, like brown dwarfs, transition between two spectral classes. An 
eventual detection of methane absorption in transmission would also help lift that 
degeneracy. 

We obtained a second series of high-resolution spectra during transit, using 
HARPS. We reanalyse the Rossiter-McLaughlin effect. The data now favour an aligned 
orbital solution and a stellar rotation nearly three times slower than stellar line broad¬ 
ening implies. A contribution to stellar line broadening, maybe macroturbulence, is 
likely to have been underestimated for cool stars, whose rotations have therefore been 
systematically overestimated. 

Key words: planetary systems - planets and satellites: atmospheres - planets and 
satellites: individual: WASP-80b - binaries: eclipsing - brown dwarfs - Hertzsprung- 
Russell and colour-magnitude diagrams 


2 Amaury H. M. J. Triaud et al. 


servatory, Chile: HARPS on the ESO 3.6m (Prog ID 089.C-0151 
& 091.C-0184), the Swiss Euler Telescope, and TRAPPIST. The 
data is publicly available at the CDS Strasbourg and on demand 
to the main author, 
f E-mail: triaud@mit.edu 


1 INTRODUCTION 

Orbital migration, disc-driven or dynamical in origin, has 
created hot Jupiters (Mayor & Queloz 1995; Lin, Boden- 
heimer & Richardson 1996; Rasio & Ford 1996). Having a 
high probability to transit, they opened up the opportu¬ 
nity to study the atmospheric content of extrasolar gas gi¬ 
ants without having to spatially resolve them (e.g. Seager & 
Deming 2010). In addition, there is hope that the chemical 
abundances of these planets carries information about where 
they formed (Oberg, Murray-Clay & Bergin 2011). This can 
be combined with those systems’ architectures (multiplicity, 
eccentricity, orbital inclination) and give a greater appreci¬ 
ation of the processes generating the diversity and width of 
parameter space that exoplanets occupy. 

Their proximity to their host star influences what we 
can know about these gas giants. For instance, water in the 
upper atmosphere of the planet will be in the gas phase as a 
result of the intense irradiation from the host star; its abun¬ 
dance can be measured. This enables us to infer the oxygen 
content, a quantity which is hard to obtain in the case of a 
cold planet like Jupiter (Madhusudhan et al. 2014). Then, 
the carbon to oxygen abundance ratio can be quantified and 
linked to the chemistry of the protoplanetary disc and ul¬ 
timately, to the formation process of the planet (Bolton & 
the Juno Science Team 2010). 

Proximity also fosters the appearance of a host of exotic 
phenomena which remain to be fully investigated. For in¬ 
stance, irradiation causes hot Jupiters to be inflated, whose 
exact cause (or causes) continues to be an active area of 
research (e.g. Demory & Seager 2011 and references there¬ 
in) . The temperatures of the majority of hot Jupiters whose 
atmospheres have been studied, are in the range 1,500 to 
2,000 K. Their densities are lower than colder planets which 
makes them practical targets for transmission spectroscopy, 
at transit. Their bloated size and elevated temperature are 
also advantages when carrying emission spectroscopy, at oc- 
cultation. 

Hot Jupiters have sizes and dayside temperatures that 
makes them resemble late M-dwarfs and L-type brown 
dwarfs. They share a common colour-magnitude space (Tri¬ 
aud 2014). These non-irradiated, self-luminous objects are 
covered by dust clouds, produced by the condensation of var¬ 
ious species such as silicates, aluminium and iron bearing 
molecules as well as titanium and vanadium oxides (Kirk¬ 
patrick 2005), which become patchy with decreasing tem¬ 
peratures (Ackerman & Marley 2001; Burgasser et al. 2002; 
Artigau et al. 2009; Radigan et al. 2012; Gillon et al. 2013; 
Buenzli et al. 2014; Crossfield et al. 2014). In the near- 
infrared bands Y, J, H, and K, the L-T transition is charac¬ 
terised by a blueward shift where the hotter, deeper parts of 
the atmosphere start becoming visible through cloud clear¬ 
ings. In the mid-IR [3.6] and [4.5 /rm] bands, that transition 
is defined by a red ward slope, produced in part by a shift of 
the blackbody peak towards cooler temperatures but also, 
by an increase in CH 4 absorption (and its associated reduc¬ 
tion in CO and CO 2 ) (Kirkpatrick 2005). 

The presence of clouds and hazes is suspected to hamper 
the detection of molecules in the transmission spectra of hot 
Jupiters (e.g. Sing et al. 2011; Jordan et al. 2013; Mandell 
et al. 2013; Pont et al. 2013). If the analogy between self 
luminous ultra-cool dwarfs and irradiated gas giants holds, 
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then planets with daysides similar to mid to late T-dwarfs 
ought to be mostly cloud free, presenting to our line of sight 
a clear terminator with strong methane absorption. In the 
two Spitzer bands we used, a T-dwarf like planet would have 
an illuminated face redder than the typical hot Jupiter when 
using the two Warm Spitzer channels (Triaud et al. 2014). 

WASP-SOb is a special planet. This gas giant, discov¬ 
ered by the WASP collaboration (Pollacco et al. 2006; Collier 
Cameron et al. 2007; Hellier et al. 2012), orbits a bright (V 
= 11.9, K = 8.4, W1 = 8.28, W2 = 8.31) late K dwarf/early 
M-dwarf. Despite a 3 day orbit, typical of hot Jupiters, it 
has an equilibrium temperature close to only 800 K (Triaud 
et al. 2013). This is one of the coolest gas giants whose at¬ 
mosphere can be studied in transmission and emission using 
current instrumentation. As such, it is an obvious target to 
observe to verify how temperature affects planetary spectra 
and to continue the exploration in colour-magnitude space. 
This is one of only a few systems expected to resemble a T- 
dwarf (Triaud et al. 2014). In addition, of the other planets 
with similar or colder temperatures that have been stud¬ 
ied, all are Neptune-like planets, orbiting M-dwarfs. WASP- 
80b fills an empty space between those small, warm Nep- 
tunes, and the more massive and large, hot Jupiters. The 
warm Neptunes mostly have featureless transmission spec¬ 
tra (Kreidberg et al. 2014; Knutson et al. 2014; Ehrenre- 
ich et al. 2014) , with the notable exception of HAT-P-llb 
(Fraine et al. 2014). WASP-80b has a larger scale-height, 
more amenable for transmission spectroscopy. 

To improve the characterisation of this fascinating 
planet, we observed some of its transits and occultations at 
3.6 and 4.5 p,m with Spitzer, and gathered new ground-based 
transit photometry and spectroscopy. We combine these new 
observations with published high-quality data, and perform 
a global analysis of the resulting extensive data set. We 
present here the results of this work. Our paper is simply 
organised: we first present our data collection and analysis, 
and then discuss our results and conclude. 


2 DATA COLLECTION AND ANALYSIS 

We used the photometric and spectroscopic data presented 
in Triaud et al. (2013) as well as the lightcurves collected 
and analysed by Mancini et al. (2014) (see Table Al and 
Fig. A3). To these, we add ground and space-based photom¬ 
etry, described in the first subsection, and new spectroscopic 
information reported in the second subsection. 

2.1 Broad band photometry at transit and at 
occultation 

2.1.1 New Sloan-z’ and Gunn-r’ transit photometry with 
TRAPPIST and Euler 

We captured three additional transits during the 2013 sea¬ 
son, using TRAPPIST (Jehin et al. 2011). One of these - on 
the night starting on 2013-06-15 - was acquired in tandem 
with the Euler telescope. They are also simultaneous with 
the DFOSC and GROND lightcurves presented by Mancini 
et al. (2014), and with the Rossiter-McLaughlin effect de¬ 
scribed in the following subsection. That night, five tele¬ 
scopes observed WASP-80b simultaneously from La Silla. 


The TRAPPIST observations made use of a Sloan z’ 
filter, with each frame exposed for 13s. The transit observed 
with Euler was acquired through a Gunn-r’ hlter, and an 
exposure time of 50s. Observations were set up similarly 
on both telescopes. Notably, we kept the positions of each 
star on the CCD chip within a box of a few pixels. This is 
achieved using a ‘software guiding’ system which regularly 
derives an astrometric solution on the science images, and 
sends pointing corrections to the mount when needed. This 
improves the photometric precision by reducing the instru¬ 
mental noise due to flat-fielding errors. Table Al presents a 
log of the observations. 

After a standard pre-reduction (bias, dark, flat-held cor¬ 
rection), stellar fluxes were extracted from our images using 
the IRAF/DAOPHOT^ aperture photometry software (Stetson 
1987). For each transit, several sets of reduction parameters 
were tested. We kept the set yielding the most precise pho¬ 
tometry on the stars having a brightness similar to WASP- 
80. After a careful selection of reference stars, differential 
photometry was then obtained. The resulting light curves 
are shown in the appendix (Fig. A3). 

2.1.2 Spitzer transit and occultation photometry at 3.6 
and 4.5 fim 

We requested time to observe WASP-80 with Spitzer (PID 
90159; Triaud) and captured six eclipses of the system with 
IRAC (Fazio et al. 2004). We obtained data covering one 
transit and two occultations at 3.6 ^m and one transit and 
two occultations at 4.5 pm. The scheduling of these observa¬ 
tions was based on the orbital solution presented by Triaud 
et al. (2013). We used the subarray mode (32 x 32 pixels, 

1.2 arcsec/pixel) with exposure times of 0.4s (3.6 pm) and 
2s (4.5 pm). The resulting sets of 64 individual subarray im¬ 
ages were calibrated by the Spitzer pipeline version S19.1.0 
and are available on the Spitzer Heritage Archive Database^. 

The data were reduced in the same way for all six time 
series. We first converted fluxes from the Spitzer units of 
specific intensity (MJy sr“^) to photon counts. We then per¬ 
formed aperture photometry on each subarray image with 
IRAF/DAOPHOT. We tested different aperture radii and back¬ 
ground annuli, and found the best results with an aperture 
radius of 2.5 pixels and a background annulus extending 
from 11 to 15.5 pixels from the centre of the point-spread 
function (PSF). We measured the centre and width of the 
PSF (FWHM) by fitting a 2D-Gaussian profile on each im¬ 
age. We then inspected at the x-y distribution of the mea¬ 
surements, and discarded the few measurements having a 
visually discrepant position relative to the bulk of the data. 
For each block of 64 subarray images, we then discarded the 
discrepant values for the measurements of flux, background, 
X and y positions, and PSF widths in the x- and p-direction, 
using a lO-cr median clipping for the six parameters. We av¬ 
eraged the remaining values, taking the errors on the average 
flux measurements as photometric errors. At this stage, we 

^ IRAF is distributed by the National Optical Astronomy Obser¬ 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

^ http://sha.ipac.caltech.edu/applications/Spitzer/SHA 
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(a) bin = 2.5 min 


(b) bin = 5 min 




Hours, from mid-transit Hours, from mid-orbit 

Figure 1. Photometric timeseries covering the transit (left panel) and occultation (right panel) of WASP-80. For visual convenience, 
the data are presented binned and vertically offset. They are labeled by central wavelength of the broad-band filter used. Un-corrected 
and un-binned data are available in the appendices. 


used a moving median filter in flux on the resulting light 
curves to discard outlier measurements due to cosmic hits, 
for example. The resulting lightcurves can be inspected in 
the appendix, in Figs. Al and A2. 

2.2 HARPS observations at transit 

On the night starting on 2013-06-15, simultaneously with 
much of the ground-based photometry, we obtained a series 
of 26 spectra with HARPS (Prog.ID 091.C-0184; Triaud) at 
a cadence of approximately 10 minutes. These new data are 
of higher quality than the transit spectroscopy presented in 
Triaud et al. (2013) as a result of better weather conditions 
(clear sky, 0.6 arcsecond seeing) and higher cadence. Seven 
additional spectra were obtained on the nights leading to 
and following the transit night. They helped further contain 
the orbital parameters. 

The data were reduced using the standard HARPS re¬ 
duction software, and the radial velocities were extracted 
by cross-correlating each spectrum with a K5 mask, as was 
done in Triaud et al. (2013). 

2.3 Global data analysis 

Strong constraints on the system parameters can be derived 
by performing a global Bayesian analysis of the photometric 
and radial velocity time series. We constructed the poste¬ 
rior probability distributions functions of the global model 
parameters using an adaptive Markov Chain Monte-Carlo 
(MCMC) code described in Gillon et al. (2014) (and refer¬ 
ences therein). 

Transits and occultations were modelled with the algo¬ 
rithm written by Mandel & Agol (2002). Each light curve is 


multiplied by a baseline model, which accounts for other as- 
trophysical and instrumental effects resulting in photomet¬ 
ric variations. A quadratic limb-darkening law was assumed 
for the transits. For each light curve, a baseline model (see 
Table Al) was selected by way of minimising the Bayesian 
Information Criterion (BIC; Schwarz 1978). The radial ve¬ 
locities were modelled by a Keplerian orbit for the star com¬ 
bined to a systemic velocity, and to the spectroscopic transit 
model of Gimenez (2006). 

For the ground-based light curves, a second-order time 
polynomial is systemically assumed for each baseline. It ac¬ 
counts for a possible low-frequency stellar variability, for a 
difference in colour between the target and its comparison 
stars, and for their consequent differential extinction. The 
modelling of other effects (dependance of fluxes on FWHM, 
on position, meridian flip) were also required for a few light 
curves (see Table Al). 

For the Spitzer photometry, our baseline models rep¬ 
resenting the well-documented ‘ramp’ and ‘phase-pixel’ ef¬ 
fects (e.g Knutson et al. 2008; Lewis et al. 2013), are similar 
to those employed in the analysis of GJ 1214 (Gillon et al. 
2014); we refer the reader to that paper for a detailed de¬ 
scription. 

The only informative prior probability distribution 
functions assumed in our analysis are for the stellar mass 
M*, effective temperature Teff, metallicity [Fe/H], and 
for the limb-darkening coefhcients of each bandpass. We 
assumed normal prior distributions ^(0.58,0.05^) Mq, 
A(4145,100^) K, and (—0.14,0.16) dex, for M*, Tes, and 
[Fe/H] respectively (Triaud et al. 2013). Normal prior dis¬ 
tributions are also assumed for the coefficients ui and U 2 
of the quadratic limb darkening law. The parameters of 
our Gaussians were interpolated from the tables of Claret 
& Bloemen (2011) for each of the bandpasses we observed 
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[h] 

Table 1. Results of our MCMC fits to the photometric and spec¬ 
troscopic data for WASP-80. Errors on the last two digits of 
each parameters, are given in brackets. Asterisks mark param¬ 
eters which are controlled by priors. 


Parameters 

Units 

Values 

the star 



Teff 

K 

4143*-+®^^ * 

[Fe/H] 

dex 

-0 13'+^®) * 

U.iO _i- 

logs* 

cgs 

4 663'+®®^ 

4.DDO,_ig. 

M* 

Mq 

0.577 +®^ * 

R* 

Rq 

0.586 +1^ 

P* 

P0 


V sin 2* 

km s~^ 

1 27'+®^) 

the planet 

P 

day 

3.06785234|+®g^^ 

To 

BID 6487.425006 +25 

D (3.6 /rm) 

- 

0.02937 +®3j 

W 

day 

0.08878 +3® 

b 

Rq 

^•^-^^(-22) 

K 

m s“^ 

109 

iuy.u^_4 4^ 

logflp 

cgs 

3 145^+®®^ 

a/Ri, 

- 

12.63+^3 

Mp 

Afjup 

0-538 -36 

Rp (3.6 /rm) 

-^Jup 

u.yyyr_3is 

Pp 

PJup 

0 539 +®® 
u.oja,_24s 

Teq 

K 

825^'*’^^^ 

a 

AU 

0.0344|+3°| 

ip 

deg 

89.02 +“ 

0 

deg 

14(+i5) 

(-14) 

e 

- 

0 002^"''^^^ 
u.uu2:(_02) 

UJ 

deg 

94(+12o) 

^ (-21) 

star to planet ratio, Rp/Ri, 


0.47 /rm 

- 

0.17006|+^' 

0.62 

- 

0.17169+®® 

0.77 Atm 

- 

0.17220 +g^ 

0.80 /rm 

- 

0 1717S'+®®1 

0.90 /rm 

- 

0.17127+®g® 

3.6 pm 

- 

0.17137^3^ 

4.5 pm 

- 

0.17230j+«5 

occultation depths 


3.6 pm 

ppm 

455(+l00) 

400(_ioo) 

4.5 pm 

ppm 

944 (+ 64 ) 

brightness temperatures 


3.6 pm 

K 

gQ]^(+68) 

4.5 pm 

K 

oqo(+58) 


in, and corresponding to a star with T^g = 4145 ± lOOK, 
logg* = 4.689 ± 0.013 dex, and [Fe/H] = -0.14 ± 0.16 dex 
(Triaud et al. 2013). No priors were applied to the Rossiter- 
McLaughlin effect^. 

A first chain of 10 000 steps was launched. It brought to 
our attention the need to rescale our photometric errors in 
order to account for an extra-white noise and for some cor¬ 
related noise. We also estimated the amount of stellar noise 
(“jitter”) in the radial velocities and added this quadrati- 
cally to the radial velocity errors (see Gillon et al. 2012 for 
details). The error rescaling factors for the light curves are 
given in Table Al. For the CORALIE radial velocities, the 
deduced jitter noise was 8.1 m s”'^; the HARPS data re¬ 
quired none. With errors corrected, we performed our main 
analysis that consisted of two chains composed of 100 000 
steps. They converged according to the statistical test of 
Gelman & Rubin (1992). Table 1 collects our derived pa¬ 
rameters and their Icr confidence region. 

The resulting light curves, corrected from systematic 
effects and and co-added according to wavelength, are dis¬ 
played in Fig. 1. The transits (Fig. la) have been binned 
in 2.5 minute bins (for visual convenience only), and the 
occultations (Fig. lb) in 5 minute bins. 

2.4 A study of the Rossiter-McLaughlin effect 

In the discovery paper, we were unable to distinguish the im¬ 
pact parameter from zero, which led the fitting procedure of 
the Rossiter-McLaughlin effect to produce degenerate solu¬ 
tions for the projected spin-orbit angle /3 (e.g. Triaud et al. 
2011; Albrecht et al. 2011). Here we study several possi¬ 
bilities, getting further into some details of the Rossiter- 
McLaughlin effect that the global analysis did not touch. 

The Rossiter-McLaughlin effect was modelled using the 
prescription of Gimenez (2006), in the same manner as in 
Lopez-Morales et al. (2014), using priors as input for all the 
parameters except those controlling the effect. Both time 
series were adjusted as a single dataset, since no offset was 
detected between them. As in Triaud et al. (2013), we in¬ 
vestigated three different propositions: 

• using no priors on the projected stellar rotation usini*; 

• using a prior, obtained from spectral line broadening, 
usinf* = 3.55 ± 0.33 km s“^ (Triaud et al. 2013); 

• applying a prior on /3 = 0 ± 10°. 

The first and third propositions produced similar, 
aligned, solutions, whereas the second converged on a in¬ 
clined orbit. Thanks notably to the Spitzer photometry, 
which estimated a non-zero impact parameter, /3 could be 
resolved (54 ±7°). However, this solution yields a poorer fit 
(Ay^ = 30.7, for the same number of degrees of freedom) 
in part because it remains in 2cr tension with the adopted 
prior, supporting incompatibility between the value mea¬ 
sured thanks to the Rossiter-McLaughlin effect and that es¬ 
timated from stellar line broadening. We thus favour the 
results from the first proposition (/I = 14 ± 14°, usini* = 
1.27 ± 0.14 km s“^) and adopt them; they can be found in 

® The Rossiter-McLaughlin effect was analysed more in detail 
separately from the global analysis. The description is provided 
in Sec. 2.4. 
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- 2-10 1 2 
Hours, from mid-transit 

Figure 2. The Rossiter-McLaughlin effect caused by WASP-80b, 
as observed by HARPS, its most likely model and its correspond¬ 
ing residuals. The radial velocities are corrected for the Doppler 
reflex motion. Additional dark dots in the centre of some data 
points indicate the new data. 



Figure 3. Simulation of the Rossiter-McLaughlin effect on a star 
similar to WASP-80. Two configurations are tested here. One 
where the star has no spots, and one where it possesses spots 
at the poles covering the upper and lower 40 degrees in latitude. 
This increases slightly the amplitude of the effect, because the 
fractional contribution to the stellar rotation hidden by the planet 
has gained in importance. 

Table 1. Our estimation of vsm^^, is thus an independent 
measurement. 

The most likely fit and the data are displayed in Fig. 2. 

3 DISCUSSION 

3.1 About the Rossiter-McLaughlin effect 

Thanks to the quality of our photometric data at transit 
and the much reduced influence of limb darkening in the 
mid-infrared, it was possible to precisely measure the impact 
parameter of WASP-80b. This lifts the degeneracy between 
usinf* and /3 noted in Triaud et al. (2013), and enables 


us to show that the planet is on a coplanar orbit instead 
of inclined, as we had previously reported. The amplitude 
of the Rossiter-McLaughlin effect implies a slower rotation 
for the star than is estimated from stellar line broadening. 
The values are incompatible at the 7cr level. In Triaud et al. 
(2013), we took this difference as a sign that the planet’s 
orbital plane was severely inclined (which the poor precision 
on the impact parameter then allowed), the data are now of 
sufficient quality to instead favour a coplanar solution. This 
brings this system more in line with the patterns of spin- 
orbit alignement with stellar host parameters proposed by 
Schlaufman (2010) and Winn et al. (2010). 

We attempted to understand this discrepancy by veri¬ 
fying how usini* from spectral line broadening is affected 
by activity, and what the resulting Rossiter-McLaughlin ef¬ 
fect would look like. To do this, we used tools developed in 
Herrero, Ribas & Jordi (2014). 

The star is chromospherically active (Triaud et al. 2013; 
Mancini et al. 2014), and we could also detect, from the 
CORALIE spectra, that the Ha equivalent width of WASP- 
80 has a wider variation epoch to epoch than intra-night 
measurements show, but presents no periodicity. Despite 
these indications of stellar activity, there is no rotational 
modulation in the photometry (Triaud et al. 2013). Some 
spot configurations can lead to weak photometric signal: 
1 ) well centred and stable polar spots, 2) a continuously 
spotted active latitude, or 3) a near homogeneously spotted 
surface. Those configurations were simulated, and produced 
slight differences in the estimation of usini*, but the dif¬ 
ferences are less than 100 m s“^. Fig. 3 shows a simulated 
Rossiter-McLaughlin effect for a spotless star, and one with 
large polar spots. The difference between both models is of 
order 1 m s“^ - well within error bars. 

Furthermore, stellar spots outside of the transit chord 
produce a slope in the transmission spectrum (Jordan et al. 
2013; McCullough et al. 2014). The important coverage nec¬ 
essary to significantly alter the Rossiter-McLaughlin effect 
would presumably leave a trace that is not detected in the 
multi-wavelength photometry presented in Fig. 5. 

Resolving this v sin discrepancy is outside of the 
range of topics this paper aimed to study, however we note 
that most laws attempting to replicate macroturbulence in 
stars suggest this phenomenon’s contribution to stellar line 
broadening decreases with stellar effective temperature, and 
vanishes for the coolest stars (e.g. Gray 2008). Doyle et al. 
(2014) have instead shown that it may plateau to a set value. 
Sadly, WASP-80 is too cool and falls out of their revised 
calibration. Our issues only highlights how tricky it is to es¬ 
timate V sin i* from stellar line broadening. It reminds us of 
the vigilance we need to keep when including these values as 
priors in an analysis of the Rossiter-McLaughlin effect. Still, 
the fact remains that WASP-80 has wide lines compared to 
others stars with a similar B — V (Triaud et al. 2013). If its 
width is equivalent to usini* = 1.3 km s“^ then all other 
late K and early M-dwarfs must be rotating much slower 
than they have been thought to do until now. As such, the 
Rossiter-McLaughlin effect may become a means to study 
the effect of macroturbulence by providing an independent 
measurement of usini*. 
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Figure 4. Two colour-magnitude diagrams highlighting the position of the dayside of WASP-SOb. Circles represent planets; GJ 436b 
and the nightside of HD 189733Ab have slightly different symbols to attract attention. Diamonds show the location of ultra-cool dwarfs 
obtained from Dupuy & Liu (2012), colour-coded as a function of spectral type. Grey symbols indicate no spectral classification. The 
blackline outlines the location of blackbodies with of a size similar to WASP-SOb, but with effective temperatures ranging from 4,000K 
to 400K. The empty diamond on this line show SOOK. WASP-SOb has one of the reddest dayside amongst all hot Jupiters, but remains 
bluer than HD lS9733Ab’s nightside. 


Table 2. Predicted magnitudes and occultation depth for various 
broad bands, and for two scenarios representing the two alterna¬ 
tives for WASP-SOb: a T4 spectrum, or a blackbody. To compute 
the magnitude for the T4 type, we used the relations produced 
in Dupuy & Liu (2012), with their errors. The uncertainty on 
the distance is propagated to provide a credible range of val¬ 
ues. Our detections fall within those boundaries. Those errors are 
only useful when considering each band in isolation. Occultations 
depth were estimated using stellar visual magnitudes obtained 
from WISE (Gutri & et al. 2013) and 2MASS (Gutri et al. 2003). 


band 

absolute 

magnitude 

visual 

magnitude 

occultation 
depth [ppm] 

J2MASS 

if a T4 

14.3 ±0.4 

17.9 ±0.5 

340 ± 180 

H2MASS 

14.1 ± 0.4 

17.7 ± 0.5 

210± 120 

K2MASS 

14.1 ± 0.4 

17.7 ± 0.5 

182 ± 98 

[3.6] 

12.9 ±0.3 

16.5 ± 0.4 

520 ±210 

[4.5] 

12.3 ±0.2 

15.9 ±0.3 

930 ± 270 

[5.8] 

12.0 ±0.3 

15.6 ±0.4 

1150 ± 470 

J2MASS 

if a 900K blackbody 

20.7 24.3 ±0.25 

1.0 ±0.2 

H2MASS 

17.5 

21.1 ± 0.25 

9.2 ± 2.2 

A' 2 MASS 

15.5 

19.1 ± 0.25 

50 ± 12 

[3.6] 

13.0 

16.6 ±0.25 

480± 110 

[4.5] 

12.3 

15.9 ±0.25 

930 ± 220 

[5.8] 

11.7 

15.3 ±0.26 

1510 ± 370 


3.2 About the colour of WASP-SOb 

Flux-drops during occultation can be combined with the 
apparent magnitudes of WASP-80, as measured by WISE 
(Wright et al. 2010) in its W1 & W2 channels (equivalent 
to Spitzer's channels 1 & 2; Triaud et al. 2014) to com¬ 
pute WASP-80b’s dayside apparent magnitudes: mp.g] = 
16.63 ±0.26 and m[ 4 . 5 ] = 15.872 ±0.078. This implies a very 
red colour for WASP-80b: [3.6]-[4.5] = 0.76±0.27, compared 
to the average colour of hot Jupiters of 0.149±0.025 (sample 
taken from Triaud et al. (2014), except for one significantly 
updated value on HD 209458b; Zellem et al. 2014). 

We reproduced the methodology applied in Triaud et al. 
(2014) and find that the WASP-80 system has a distance 
modulus m — M = 3.55 ± 0.26 (51.3 ± 6.2 pc). The absolute 
magnitudes of WASP-80b’s dayside are therefore: Mp gj = 
13.08±0.37 and M[ 4 , 5 ] = 12.32±0.27. This is slightly fainter 
than the nightside of HD 189733Ab; it is also bluer in those 
specific bands (Knutson et al. 2012; Triaud 2014). 

The Spitzer occultations correspond to brightness tem¬ 
peratures of 9011®2 888^^ K for channel 1 and 2 respec¬ 

tively, slightly above its equilibrium temperature of SOOK. 
These values have been computed using the stellar models, 
computed by Kurucz (1993), that were compatible with the 
stellar parameters that we have used for WASP-80. The er¬ 
ror budget contains the observational uncertainties on the 
occultations’ depths and on the stellar parameters, as in De- 
mory et al. (2013). 

Dupuy & Liu (2012) present relations between abso¬ 
lute magnitudes in certain bandpasses, including those we 
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observed in, and spectral classification of ultra-cool dwarfs. 
The dayside of WASP-80b corresponds to a T4 dwarf, in 
both the [3.6] and [4.5] channels. T4 is the latest type con¬ 
sidered within the L-T transition. The T4 type comprises 
effective temperatures ranging between 1100 and 1300 K 
(Stephens et al. 2009). A 200 K increase can be compen¬ 
sated by a 30% reduction in radius and produce a similar 
luminosity. According to Baraffe et al. (2003), brown dwarfs 
can be 0.75 i?jup, which would make WASP-80b’s dayside 
fit with a T4’s flux, within uncertainties. Provided that this 
is indeed the case, our results may imply that for planets, 
the L-T transition happens at cooler temperatures than for 
brown dwarfs^. 

We do not intend here to attach a spectral classification 
to WASP-80b, but only to provide a reference point. These 
relations have been built from -and are therefore only valid 
for- ultra-cool dwarfs: self-luminous objects, whose stronger 
gravities impact the structure and type of clouds layers. 
While for brown dwarfs, a change in composition leads to 
small variations in the emerging spectrum, planetary spectra 
can be severely affected. This is particularly relevant since 
planets cover an important range in metal enrichment as 
seen in the Solar system. For similar composition, irradia¬ 
tion changes a planet’s spectrum dramatically too (Burrows 
2014) and different chemistry can be expected between the 
night and day sides of a planet, rendering the classification 
of a planet arduous at best (Cooper & Showman 2006). In¬ 
deed, the phase curve of HD 189733Ab presents evidence 
that its atmosphere goes through a change in colour similar 
to a spatial transition between L to T (Triaud 2014). 

Our values were placed into two colour-magnitude dia¬ 
grams (Fig. 4), presenting the same data as in Triaud et al. 
(2014) (except for HD 209458b; Zellem et al. 2014 and the 
nightside of HD 189733Ab; Triaud 2014). WASP-80b sits 
just on the T-dwarf sequence, which confirms our earlier 
statement about its resemblance to a T4 spectral type. The 
redward direction of the T-dwarf sequence (in the [3.6]-[4.5] 
colour) is characterised by the appearance of strengthening 
CH 4 absorption in the 3.6 pm channel ( and a correspond¬ 
ing decrease in CO: see Patten et al. 2006; Zahnle & Marley 
2014). 

Two blackbody loci, corresponding to a given size of 
1 7?jup, for all effective temperatures between 4,000 and 
400 K, are also drawn on Fig. 4, with a mark where 800 K 
is. WASP-80b coincides with blackbodies of temperatures 
890lggK and 900t|oK for [3.6] and [4.5 pm] respectively. 
These temperatures are equivalent to brightness tempera¬ 
tures but are estimated differently than the values provided 
in Table 1. In this new approach, we make no assumption 
on the shape of the stellar photosphere, but instead are de¬ 
pendent on the assumed distance. It is reassuring to see 
agreement between these two different methods.. 

In the [3.6]-[4.5] colour, the T-dwarf sequence and the 
blackbody both slope toward redder colours. This is in con¬ 
trast to the near IR J, H Ks and mid-IR [5.8 pm] 
bands where the L-T transition is characterised by a blue- 
ward slope (while the blackbody remains redward; Triaud 


^ that transition is expected to occur around 500K (Zahnle &; 
Marley 2014). However their work was carried for self-luminous 
objects only. 


et al. (2014)). Predictions are presented in Table 2. Observ¬ 
ing in the J, H and Ka bands is particularly promising to 
distinguish between both propositions. 

The brown dwarf sequences and the blackbody expec¬ 
tations are furthest apart, over the two Warm Spiteer chan¬ 
nels, for planets with equilibrium temperatures intermediate 
between WASP-80b and 1,000K (e.g. WASP- 8 b; Cubillos 
et al. 2013). Observations of additional systems in this range 
will help interpret the results for WASP-80b. Those obser¬ 
vations would also verify the existence of diversifying atmo¬ 
spheric properties with increasing magnitude, as revealed by 
an increase in the range of colour occupied by cooler planets 
(Triaud et al. 2014). 


3.3 About the transmission spectrum 

The planet-star radius ratio (i?p/i?*) in each of our broad 
bandpasses can be found at the bottom of Table 1. They 
are graphically presented in Fig. 5, and except for the 
value at 0.47 pm, are all consistent with each other. The 
inverse-variance weighted average of these values is i?p/7Z* = 
0.17167 ± 0.00020, which is also plotted on Fig. 5. The data 
agree well with a horizontal line: = 8.3 ±4.1, which for 

six degrees of freedom (seven wavelengths and one free pa¬ 
rameter) makes a reduced x? = 1.4 ± 0.7. We see no rea¬ 
son to compare our data to more complex models (even 
a slope would be too much). Our analysis refines the re¬ 
sults of Mancini et al. (2014). WASP-80b joins a number 
of hot Jupiters whose spectra are approximately fiat within 
measurement errors (e.g. Pont et al. 2013; Mancini et al. 
2013). Flat transmission spectra have also been observed for 
warm Neptunes such as GJ 1214b (Berta et al. 2012; Kreid- 
berg et al. 2014), GJ436b (Knutson et al. 2014), GJ3470b 
(Ehrenreich et al. 2014), a result that has been interpreted 
as indicating the presence of clouds or as a sign of high 
metal content. 

Fukui et al. (2014) also observed WASP-80b. They find 
Rp/Ri, = 0.17000 ± 0.00040 (averaged over the J, H and 
Ka bands, who are all consistent with each other) a value 
similar to the measurement at 0.45 fim (Sloan g’) shown in 
Fig. 5. In the visible they measure Rp/Ri, = 0.1741 ±0.0010, 
which is higher than what we find®. In their analysis, some 
transit parameters were kept fixed which usually results in 
underestimated error bars and would make any dispersion 
appear statistically more significant than it really is. For this 
reason, we have not included their results in our analysis. 

If the redness of WASP-80b’s dayside is an indication 
of methane absorption then a corresponding signature is to 
be expected in transmission. Fortney et al. (2010) present 
models for planets with equilibrium temperatures of 1000 K 
and 500 K. While the former remains mostly influenced by 
H 2 O, Na and K absorption, the latter shows a transition to 
an atmosphere dominated by CH 4 absorption, closer to what 
we expect in the case of WASP-80b. This simple picture can 
be considerably altered in the presence of an opaque cloud 
cover (possibly assembled from soots, created through via 
photochemical reactions with methane (Moses et al. 2011)). 


® those values only account for their transit depths, after remov¬ 
ing those of Mancini et al. (2014) which they had included too. 
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Figure 5. Ratios of the planet to star radii as a function of wavelength. The dotted line is the weighed mean, and the grey area the error 
on that mean. A model atmosphere (H 2 : 93% , He: 7%, H 2 O: 10“^ , CH 4 : 10“^ , CO: 10“'^ , Na: 10“® and K: 10~^) is also displayed, with 
two different resolutions, drawn using different shades. This model is not fitted, and is here purely for illustrative purposes. Coloured bars 
correspond to the measurements as expected from the model; their widths illustrate the full width at half maximum of their respective 
filters. 


For illustrative purposes, we calculated a model spec¬ 
trum for WASP-SOb, which we display in Fig. 5 (Ehrenre- 
ich et al. 2006, 2014)®. It was set with the following abun- 
dances: Hz: 93% , He: 7%, H 2 O: 10-^ CH 4 : 10“^, CO: 
10“'^, Na: 10“® and K: 10“^ as in Moses et al. (2011). 
The visible regions are dominated by alkali absorption, with 
some Rayleigh scattering, produced purely from Hz. We note 
that the bluest measurement is particularly discrepant to 
that model. A visual comparison of that model to our data 
also shows that one of the strongest methane features peaks 
around 3.5 ^m. Methane should enhance a planet’s radius at 
[3.6] compared to [4.5], which we do not see in transmission 
( although we detect an equivalent effect in emission). This 
is the most we can reliably affirm. 


4 CONCLUSION 

WASP-SOb is on an aligned orbit and is the second faintest 
planet with an emission measured in each of the Warm 
Spitzer channels, after GJ 436b at [3.6] (Lanotte et al. 2014) 
and 55 Cnc e at [4.5] (Demory et al. 2012), which are indi¬ 
cated in Fig. 4. It is the faintest measured in both bands and 
for which a colour can be computed. The planet’s dayside 
is amongst the three reddest daysides of any hot Jupiter. 
HD 189733Ab’s nightside remains the measurement furthest 
to the red (Knutson et al. 2012; Triaud 2014). 

GJ 436b and WASP-SOb have similar equilibrium tem¬ 
peratures (Stevenson et al. 2010) and both show a flat spec¬ 
trum (Knutson et al. 2014). However, we find that their 

® nota bene: this model does not solve to radiative equilibrium 


colours are radically different which can be interpreted as 
evidence for different chemical compositions (as proposed 
in Stevenson et al. (2010); Lanotte et al. (2014)). We note 
here that GJ 436b has a radius significantly smaller than 
WASP-SOb’s; one could wonder if both ought to be com¬ 
pared since gravity has an important role to play on the 
condensation of chemical species into clouds layers (GJ 436b 
is 2.3 times denser than WASP-SOb). However, we find that 
there is some credence to that: a change in radius will cause 
a change in absolute magnitude. Correcting GJ 436b’s ra¬ 
dius to WASP-SOb’s (1 i?jup) amounts to comparing abso¬ 
lute magnitude produced by similar areas. If it were that 
large, GJ 436b would be 2.1 magnitude brighter. Using the 
values presented in Triaud et al. (2014), this would mean a 
Af[ 3 , 6 ] = 13.13 ± 0.28, on a par with WASP-SOb. Both plan¬ 
ets emit a similar flux density at 3.6 fim (while remaining 
very different at 4.5 /rm). 

WASP-SOb’s atmosphere is ambiguous. While WASP- 
SOb presents a dayside with a luminosity and a colour rem¬ 
iniscent of a mid T-dwarf (who have a patchy cloud cover 
and methane absorption), its transmission spectrum appears 
featureless (implying overcast, without direct evidence for 
methane). The planet is also entirely compatible with a 
blackbody. Upfront the current data would fit well with 
WASP-SOb being a dusty isothermally radiating ball. How¬ 
ever, some brown dwarfs too are compatible with black- 
body colours and would too present a flat transmission. In¬ 
stead, direct spectroscopy shows they are far from resem¬ 
bling blackbodies. We need to keep entertaining the possi¬ 
bility that neither is WASP-SOb. 

Is WASP-SOb a T4? Likely not. If we have learned 
anything about exoplanets, is that they beat our simplest 
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explanations. However its current position in the colour- 
magnitude diagram permits us to raise interesting questions. 
These diagrams become a reading grid. More observations 
are dearly needed to hrst, verify WASP-80b’s colour, espe¬ 
cially in the 3.6 ^m band where our detection is only at 
the 4(7 level and secondly, to scrutinise the planet’s size as 
a function of wavelength in the near infrared, in search for 
some weak methane absorption. Alternatively, occultation 
measurements in the near infrared would confirm whether 
the planet obeys a T dwarf behaviour, constraint its flux to 
a blackbody, or uncover its own planetary behaviour. 

The methods developed by Snellen et al. (2010) to de¬ 
tect CO and H 2 O in absorption and in emission in exoplan¬ 
ets (e.g. Rodler, Lopez-Morales & Ribas 2012; Birkby et al. 
2013) could in principle be used to probe for the presence 
of methane. Observations have been attempted at the VLT 
using CRIRES targeting a region containing a forest of CH 4 
lines, around 2.3 fim. However bad weather and schedul¬ 
ing issues meant that fewer transits were obtained than re¬ 
quested. CRIRES has now been un-mounted from the tele¬ 
scope for upgrade. 

Would a transmission spectrum turn out to show no 
spectral signatures, we could postulate the following: 

• WASP-80b is a featureless blackbody, presumably be¬ 
cause of a dusty/cloudy atmosphere, that just happened to 
coincide with the T-dwarf sequence. 

• WASP-80b is part of a diverse population of warm plan¬ 
ets enveloped by a greater variety of atmospheres than 
brown dwarfs have, as proposed in Triaud et al. (2014). 
WASP-80b’s position at the conjunction between blackbody 
and T-dwarfs would therefore be a coincidence. If indeed 
the case, it implies that its composition cannot be inferred 
from comparing with ultra-cool dwarfs, such as T-dwarfs, 
whose atmospheric chemistry and optical properties would 
share little in common with planets (e.g. Zahnle & Marley 
(2014)), especially the irradiated kind. 

• WASP-80b’s colour is an integrated view of its dayside 
surface which may be in-homogenous (like Kepler-7b; De- 
mory et al. 2013). The hottest parts are possibly covered in 
clouds (like L dwarfs and early type T-dwarfs), while the 
colder zones of the dayside may remain free, clouds having 
rained out. Equatorial wind could move enough of the clouds 
to the night side, traversing the terminator and affecting the 
resulting transmitted spectrum. 

To complicate further, the range of colour indices cov¬ 
ered by the colder gas giants may be produced by clouds cov¬ 
ers having distinct optical properties, for instance, if trace 
elements specific to some planet react with the surfaces of 
suspended cloud particles, leading to different albedos and 
opacities. In addition, extra diversity can be caused by plan¬ 
ets possessing a different number and alternation of belts and 
zones with contrasting albedos (Jupiter and Saturn are dis¬ 
similar) . 

Most of these considerations rely on the fact that brown 
dwarfs and planets can be compared to each others directly. 
The exercise is currently useful since the brown dwarfs are 
well understood and that the planets are not. Eventually, 
we will come to study each population in isolation once 
planetary-specific behaviour has been identified. Reasons 
range from the irradiated conditions of hot Jupiters, to be¬ 
ing in a different gravity range compared to brown dwarfs. 


Chemical transitions very likely differ between brown dwarfs 
and planets. Gravity influences the chemical balance be¬ 
tween CO and CH 4 (Zahnle & Marley 2014). This would 
blunt the sharp L-T transition and cause a shallower slope 
in the [3.6]-[4.5] colour than for brown dwarfs, making it 
closer to a blackbody. 

Currently, it remains unclear whether WASP-80b’s at¬ 
mosphere contains any methane. Few planets have been ob¬ 
served near the L/T transition. Apart from WASP-80b, most 
are young, directly imaged systems which despite distinctive 
physical parameters, present similar atmospheric character¬ 
istics: they are close to blackbodies (in near IR; Triaud et al. 
(2014)), likely cloudy (Marley et al. 2012), and offer methane 
absorption lower than expected (Konopacky et al. 2013). 
As such, deepening the study of WASP-80b will prove im¬ 
portant. Parallel properties between the massive, far flung 
giants orbiting a young A star, and the light, mildly irra¬ 
diated giant orbiting an M dwarf would inform us about 
widely occurring phenomena in exoplanets (be they related 
to atmospheric physics or to formation processes). 

Doubtless to say, other systems intermediate between 
the general population of hot Jupiters and WASP-80b would 
be fascinating to observe, in order to study where they land 
in colour-magnitude diagrams. Observing the shape of their 
transmitted spectra can bring to the fore a relation between 
colour and chemical composition/clouds presence. 

NOTA BENE 

Dates are given in the BJD-UTC standard. The radii we 
used for Jupiter and the Sun are the volumetric mean radii. 
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APPENDIX A: 

Here are located the graphical representation of each pho¬ 
tometric time series that has been used in onr analysis, in¬ 
cluding the corrections given to the data. Table Al contains 
the journal of observations and type of corrections applied 
to the data. 
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Date 

Instrument 

Filter 


Np 

Baseline 

function 

CF 

Eclipse 

type 

7 May 2012 

TRAPPIST 

Sloan z’ 

lOs 

728 

p(t^) + 0 

2.08 

transit 

26 Jul 2012 

Euler 

Gunn r’ 

60s 

208 

+ f'whm?') 

1.36 

transit 

10 sep 2012 

TRAPPIST 

Sloan z’ 

10s 

705 

p(t^) 

1.15 

transit 

13 Jun 2013 

TRAPPIST 

Sloan 2 ’ 

13s 

749 

p(t^) + 0 

1.20 

transit 

16 Jun 2013 

TRAPPIST 

Sloan 2 ’ 

13s 

782 

v(p + + 0 

1.49 

transit 

16 Jun 2013 

Euler 

Gunn r’ 

50s 

175 

p(t^) 

2.86 

transit 

16 Jun 2013 

2.2m/GROND 

Sloan 2 ’ 

60s 

157 

p(t^) 

2.64 

transit 

16 Jun 2013 

Danish/T)FOSC 

Bessell I 

60s 

200 

p(t^) 

1.05 

transit 

16 Jun 2013 

2.2m/GROND 

Sloan r’ 

60s 

156 

p(t^) 

2.14 

transit 

16 Jun 2013 

2.2m/GROND 

Sloan i’ 

60s 

156 

p(t^) 

1.05 

transit 

16 Jun 2013 

2.2m/GROND 

Sloan g’ 

60s 

162 

p(t^) 

0.76 

transit 

3 Jul 2013 

Spitzer 

3.6 pm 


548 

p(xy^ + (1) + BM 

1.35 

occultation 

13 Jul 2013 

Spitzer 

4.5 pm 


123 

p{fwhm^ -|- P) 

1.06 

transit 

18 Jul 2013 

Spitzer 

4.5 pm 


121 

pixy^ + l^) 

0.68 

occultation 

24 Jul 2013 

Spitzer 

3.6 pm 


569 

p{fwhml, + xy^ + P) + BM 

1.78 

occultation 

27 Jul 2013 

Spitzer 

4.5 pm 


123 

p(fwhm^ + xy"^ A- P) 

1.14 

occultation 

1 Aug 2013 

Spitzer 

3.6 pm 


579 

p(fwhm^ + xy"^ + P) -1- BM 

1.33 

transit 

1 Aug 2013 

TRAPPIST 

Sloan 2 ’ 

13s 

512 


1.23 

transit 


Table Al. WASP-SOb photometric eclipse time-series used in this work. For each light curve, this table shows the date of acquisition, 
the used instrument and filter, the exposure time Texp, the number of data points, the baseline function selected for our global analysis 
(see Sec. 2.3), the error correction factor CF used in our global analysis, and the nature of the eclipse (transit or occultation) For the 
baseline function, p(e^) denotes, respectively, a A-order polynomial function of time (e = t), the logarithm of time (e = /), x and y 
positions (e = xy)^ full-width at half-maximum FWHM (e = fwhm) and FWHM in the a:-direction (e = fwhrrix)', o denotes an offset 
at the time of a meridian flip of TRAPPIST (see Gillon et al. 2012); BM denotes the use of the BLISS-mapping technique (Stevenson 
et al. 2012) to improve the modeling of the position effects (see Gillon et al. 2014). 
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Hours, from mid-orbit Hours, from mid-orbit 




Hours, from mid-orbit Hours, from mid-orbit 

Figure Al. Photometry at occultation, as extracted from the Warm Spitzer frames, in both channels, with the initial ramp removed. 
The full model to the data is shown as a plain line; the residuals are displayed underneath. The corrections, including the variations 
in intra-pixel sensitivity, are isolated and drawn as a dotted line (barely noticeable here due to the weakness of the occupation signal 
compared to the corrections). 




Hours, from mid-transit Hours, from mid-transit 


Figure A2. Photometry at transit, as extracted from the Warm Spitzer frames, in both channels, with the initial ramp removed. The 
full model to the data is shown as a plain line; the residuals are displayed underneath. The corrections, including the variations in 
intra-pixel sensitivity, are isolated and drawn as a dotted line. 
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Hours, from mid-transit 


Hours, from mid-transit 


Figure A3. Photometry at transit, in the visible wavelengths. The full model to the data is shown as a plain line; the residuals are 
displayed underneath. The corrections, for instance airmass, or changes in seeing or pointing, are isolated and drawn as a dotted line. 


















